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RELAXATION  BEHIND  STRONG  SHOCK  WAVES  IN  AIR 

This  report  summarizes  work  performed  during  a two  year  National 
Research  Council  Postdoctoral  Research  Associateship  appointment 
which  ended  30  June  1974. 

This  task  was  initiated  to  provide  a flow  field  model  for  the 
relaxation  region  of  an  explosion  shock  wave  in  air  to  be  used  in 
an  investigation  of  airshock  radiant  light  output;  in  particular, 
the  model  was  developed  for  use  in  studying  the  appearance  of 
non-luminous  airshocks  produced  by  chemical  explosions  observed 
by  Dr.  L.  Rudlin,  the  scientific  advisor  for  this  project. 

The  task  was  divided  into  two  steps.  The  first  step,  the  development 
of  the  airshock  relaxation  region  model,  was  accomplished;  however, 
time  did  not  permit  significant  work  on  the  second  step  of  the  task, 
the  formulation  of  the  radiation  model. 
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RELAXATION  BEHIND  STRONG  SHOCK  WAVES  IN  AIR 
1 . INTRODUCTION 

Physical  and  chemical  relaxation  phenomena  associated  with  a strong 
shock  wave  in  air  are  discussed  in  this  report.  The  investigation 
was  performed  for  future  evaluation  of  mechanisms  whereby  radiant 
light  is  produced  in  an  airshock.  RUDLIN  (1970)*  observed  nonluminous 
airshocks  produced  by  chemical  explosions  (HE)  for  Mach  numbers  on  the 
order  of  20  at  sea  level  conditions.  This  absence  of  expected  radia- 
tion corresponding  to  post-shock  equilibrium  conditions  prompted  this 
investigation . 

The  work  discussed  here  is  concerned  with  flow  relaxation  effects 
exclusive  of  radiation  considerations . The  philosophy  is  to  provide 
a computational  model  for  the  flow  field  which  may  then  be  used  to 
determine  the  radiation  field  separately.  This  technical  report 
compares  calculations  for  shocked  equilibrium  conditions,  ionizational 
relaxation  times,  C>2  and  N2  vibrational  relaxation  times,  and  O2 
half-dissociation  lengths  with  experimental  data.  Computations  for 
the  following  flow  conditions  are  discussed:  T = 300 °K,  P = 5 Torr, 

and  a Mach  number  range  M = 7-36 . The  pressure  of  5 Torr  is  selected 
so  that  calculations  can  be  compared  with  available  experimental  data. 
The  Mach  number  range  is  chosen  to  conform  with  the  flow  conditions 
investigated  by  RUDLIN  (1970) . 

The  airshock  chemistry  model  and  the  computer  code  used  for  this 
study  were  originally  developed  by  SHERMAN  (1968)**.  The  present  air- 
shock relaxation  model  is  introduced  by  listing  the  important 
assumptions . 

(1)  One-dimensional  steady  flow  assumed. 

(2)  Effects  of  viscosity,  thermal  conductivity,  and  molecular 
transport  neglected. 


^References  are  listed  in  7.  REFERENCES  of  text  in  alphabetical 
order. 

**The  author  wishes  to  acknowledge  several  helpful  conversations  with 
Dr.  Sherman  concerning  his  model  and  computer  code. 
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(3)  Effects  of  radiation  on  flow  field  neglected. 

(4)  Translational  equilibrium  of  all  chemical  species  assumed, 
exception  - separate  free-electron/bound-electron  equilibrium  allowed. 

(5)  Internal  equilibrium  assumed  for  rotation  and  vibration; 
exceptions  - separate  vibrational  relaxation  allowed  for  O2  and  N2 . 

(6)  Specific  thermodynamic  and  chemical  model  selected. 

The  Sherman  airshock  chemistry  model  includes  bremsstrahlung  and 
free-bound  radiative  recombination  of  N+  and  0+.  The  radiation  model 
is  not  considered  here . 

Improvements  incorporated  into  the  airshock  relaxation  model 
include  vibrational  relaxation  of  C>2  and  N2  and  separate  electron/atom 
temperature  determinations.  Also,  the  Sherman  airshock  chemistry 
model  (1968)  has  been  updated  with  more  recent  (1972)  thermodynamic 
and  chemistry  data. 

LIN  and  TEARE  (1963)  present  a theory  for  ionizational  relaxation 
which  agrees  with  experimental  data  for  the  Mach  region  M = 7-20*. 
BIBERMAN  and  YAKUBOV  (1965)  provide  a theory  for  M>29  and  ZHELEZNYAK 
and  MNATSAKANYAN  (1968)  consider  the  region  M = 26-36  as  do 
ZHELEZNYAK,  MNATSAKANYAN,  and  YAKUBOV  (1971) . They  each  cover  the 
important  transition  region  M = 26-29  within  which  the  primary  mech- 
anism for  determining  ionizational  relaxation  changes  from  control  by 
collisions  of  heavy  particles  to  control  by  electron-heavy  particle 
collisions.  The  present  theory  is  applicable  for  the  Mach  number 
regime  M = 7-36. 

Results  from  the  present  calculations  are  compared  with  available 
experimental  data.  The  results  are  summarized  as  follows: 

1.  Excellent  agreement  obtained  for  post-shock  equilibrium 
conditions,  M = 7-36. 

2.  Ionizational  relaxation  times  fall  well  within  the  scatter  of 
data  for  M = 7-20.  However,  there  is  only  fair  agreement  for  the 
range  M = 26-36. 

3.  O2  dissociation  lengths  follow  available  data  for  M = 8-16. 

4.  Vibrational  relaxation  times  for  O2  and  N2  are  reasonable  for 
the  range  M = 7-26.  However,  the  model  is  inadequate  for  N2 
vibrational  relaxation  for  M > 26. 


* The  theory  of  LIN  and  TEARE  (1963)  is  valid  to  M v 26  with  respect 
to  ionizational  relaxation  since  it  is  believed  that  the  same 
ionization  mechanism  is  dominant.  However,  there  is  no  experimental 
data  for  comparison  between  M = 20-26 . 
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2 .  THEORY 

Air  is  assumed  to  contain  11  species,  02 , N2 , 0,  N,  NO,  02+,  N2+, 
0+,  N+,  N0+  and  e.  The  pre-shock  species  are  N2(79%)  and  02(21%). 
Applying  conservation  equations  for  charge  neutrality  and  oxyaen  and 
nitrogen  atoms  reduces  the  number  of  independent  species  conservation 
equations  from  11  to  8.  The  eight  equations  chosen  may  be  expressed  in 
the  form* 


_d 

dx 


(X.) 


-1 

2nu 


(\jXKj) 


i,  j 1,...,8 


where  represents  the  mole  fraction  of  the  chemical  species  i and 
Ki  is  the  rate  of  production  of  particles  of  species  i per  unit  volume. 
Matrix  A^j  and  the  50  reaction  rates  used  are  presented  in  APPFNDIX  A. 

The  production  terms  are  defined  in  the  manner  described  by 
SHERMAN  and  KIRK  (19  65)  . 

The  types  of  chemical  reactions  considered  are: 

1.  Dissociation  - association 

2.  Rearrangement 

3.  Ionization  - recombination 

4.  Positive- Ion  charge  exchange 


The  following  reaction  rates  are  evaluated  at  the  local  translational 
(atom)  temperature:  dissociation,  association,  rearrangement, 

neutral-neutral  collisional  ionization,  and  positive-ion  charge 
exchange.  The  reaction  rates  for  electron-neutral  collisional  ioni- 
zation and  electron-ion  recombination  are  evaluated  at  the  local 
electron  temperature.  The  dissociation  rates  for  N2  and  O2  are 
corrected  for  vibration-dissociation  coupling  usinq  factors 
developed  by  HAMMERLING,  TEARE,  and  KIVEL  (1959).  See  APPFNDIX  A. 

Vibrational  relaxation  is  considered  for  species  02  and  N2 . The 
relaxation  equation  for  vibrational  energy  for  each  species  contains 
the  following  contributions. 


*See  6.  LIST  OF  SYMBOLS  of  text  for  complete  definition  of  terms. 
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dE 

v 


= Collisional  - Vibrational  energy  exchange 


+ 


Dissociational  - Vibrational  energy  exchange 


Electron  - Vibrational  energy  exchange  (N2  only) 


Collisional-vibrational  energy  exchange  considers  the  direct 
conversion  of  the  energy  of  the  colliding  particles  into  excitation 
of  the  vibrational  energy  modes.  The  characteristic  relaxation  times 
required  for  computing  this  term  for  N2  and  O2  are  taken  from  empirical 
fits  to  experimental  data  provided  by  MILLIKAN  and  WHITE  (1963) . The 
second  term  couples  the  effects  of  the  dissociation  process  to  the 
vibrational  relaxation  rate.  Vibrational  energy  is  lost  or  gained  for 
each  dissociation  or  association  event,  and  the  number  of  oscillators 
changes  during  the  dissociation-association  process  (TREANOR  and 
MARRONE  (1962) ) . The  third  term  is  only  important  for  N2  and  represents 
energy  exchange  between  free  electrons  and  the  molecular  oscillators. 
This  term  will  be  discussed  in  connection  with  the  electron  energy 
equation.  See  APPENDIX  B for  the  explicit  formulation  of  this 
vibrational  relaxation  equation. 

Vibrational  energy  exchange  between  N2  and  O2  species  is  not 
considered  in  this  model.  GENERALOV,  LOSEV,  and  OSIPOV  (1964)  showed 
that  the  major  effects  of  vibrational  energy  transfer  occur  for 
moderate  Mach  numbers  corresponding  to  final  equilibrium  temperatures 
less  than  3000°K.  TURCHAK  (1968)  considers  a Mach  number  of  M = 12 
for  an  ambient  pressure  of  10~3  atmospheres  to  be  an  upper  limit  for 
needing  to  include  vibrational  energy  exchange.  However,  Turchak  (as 
do  Generalov,  Losev,  and  Osipov)  neglects  dissociation  and  ionization 
processes.  Since  the  vibrational  relaxation  time  for  N2  is  of  the 
same  order  of  magnitude  for  these  flow  conditions  as  the  dissociation 
time  for  O2,  the  above  stated  conclusions  represent  overestimates  for 
the  effects  of  vibrational  energy  exchanges.  For  this  reason 
vibrational  energy  exchange  between  N2  and  O2  is  neglected  for  the 
present  calculations  which  extend  down  to  M = 7 . 

The  electron  energy  equation  is  given  as 


dT 


T dn 


e 1 


e e 


dx 
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which  is  a modified  version  of  the  general  expression  presented  by 
BIBERMAN  AND  YAKUBOV  (1965)  . The  differences  include  substitution  of 
heavy  particle  velocity  for  electron  velocity  and  deletion  of 
negligible  terms  accounting  for  convection  and  electron  thermal 
conduction  (BIBERMAN,  MNATSAKANYAN , and  YAKUBOV  (1971)). 


The  latter  investigators  are  able  to  reduce  the  electron  energy 
equation  still  further  to  the  "quasistationary  form"  in  which  the 

dTe  dne  ' 

terms  containing  the  derivatives  and  are  nealected.  What  is 


5 

left  is  the  term  which  represents  the  local  electron  energy 

1 


balance.  The  quasistationary  form  of  the  equation  is  not  appropriate 
for  the  present  calculation;  the  derivative  terms  are  not  negligible. 


Five  energy  flux  terms  are  included  in  the  electron  energy  balance: 

- Energy  gained  by  electrons  through  elastic  collisions  with 
ions 

Q2  - Energy  lost  by  electrons  through  electron-neutral  collisional 
ionization  reactions 

Q3  - Energy  gained  by  electrons  through  inelastic  collisions  with 
N2  molecular  oscillators 

Q4  - Energy  gained  by  electrons  from  associative  ionization 
reactions 

Q5  - Energy  lost  by  electrons  maintaining  the  quasistationary 
population  of  excited  states. 

APPENDIX  C provides  more  detailed  information  on  the  electron  energy 
equation . 


In  addition  to  the  above  set  of  equations , the  integrated  equations 
for  conservation  of  mass,  momentum,  and  energy  must  be  satisfied  and 
an  equation  of  state  is  required.  The  conservation  equations  written 
in  a coordinate  system  moving  with  the  shock  front  at  velocity  U are 
written  as 

pu  = pQ  U 

p + pu2  = pQ  + Pq  U2 

„2 

h (T , TV02 , TVN2,  T , X ) + — = h(T  , 

e i 2 o 


(X.)  ) + 2_ 

10  2 
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Enthalpy  is  defined  in  APPENDIX  D.  The  equation  of  state  for  a 
partially  ionized  gas  mixture  with  a separate  electron  temperature 
is 

p = n k (T  + Xe  (Te  - T)  ) . 
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3 . NUMERICAL  METHOD 


Conventional  numerical  methods  such  as  Runge-Kutta  or 
predictor-corrector  are  not  adequate  for  integrating  the  ordinary 
differential  equations  of  one-dimensional  chemical  nonequillibrium 
flow.  The  maximum  allowable  step  size  needed  for  such  integration 
using  these  techniques  is  very  small  and  extensive  machine  time  can 
be  required  for  the  calculations . The  step-size  requirement  of  the 
conventional  methods  is  due  to  the  "stiff"  nature  of  the  differential 
equations.  This  term  is  defined  as  follows  (SHERMAN  (1968));  the  i^ 
equation  of  the  system 

dyi 

3x“  = fi(yl'  y2 ' * * *yn'  X)  i=l , . . . ,n 

is  said  to  be  stiff  in  yj  if  a small  change  in  yj  causes  a large 
change  in  f£. 

Several  methods  for  handling  stiff  equations  with  discussions  of 
step  size  control  and  stability  criteria  are  described  in 
SHERMAN  (1968)  and  LOMAX  and  BAILEY  (1967) . The  technique  used  to 
solve  the  eight  species  production  equations  simultaneously  is  the 
two-term  complete  matrix  method  with  local  linearization  applied 
which  is  now  described  briefly. 

The  original  differential  equation,  shown  above,  is  replaced  by 
the  simple  approximate  form*  (SHERMAN  (1968)) 


dy: 

dx 


- -L, 


r^pi j (Yj  " Yix)  + Ai  i,j=l,...,n 


with  the  following  definitions 


point  1.  x = xx  , y • = y 


jl 


point  2.  x = x-l  + h , y^2  = Yji  + hfji 


dy  • 


f . , = -t— ^ - at  x = x.  , y . = y , y = y . , j 
lk  dx  k J j Jm  Jmi' 


* The  approach  to  local  equilibrium  via  exponential  decay 
dy/dx  = -P(y-yeq)  suggests  the  approximate  form. 
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dY  i 

Ai  “ at  x ■ v yd 


yii 


In  determining  Pi j , the  column  is  determined  by  permitting 

only  yj  to  vary  and  holding  all  other  y^  (m^j)  to  the  initial  values. 

After  changing  to  the  simple  approximate  form,  the  equations  are 
then  transformed  into  canonical  form  by  diagonalizina  the  matrix  P^_: 
assuming  distinct  eigenvalues.  The  transformed  system  of  equations^ 
is  solved  and  then  the  solution  in  terms  of  the  original  variables  is 
recovered  using  the  inverse  transformation.  The  matrix  methods  used 
to  accomplish  this  are  described  in  SHERMAN  (1968)  along  with  a 
discussion  of  how  complex  eigenvalues  and  identical  eigenvalues  are 
handled. 


The  differential  equations  for  vibrational  relaxation  of  O2  and 
N2  and  electron  energy  are  solved  separately  using  the  same  approxi- 
mate form.  Matrix  methods  are  not  required  since  the  three  equations 
are  uncoupled  after  they  are  linearized.  However,  the  order  in  which 
the  equations  are  solved  should  be  noted.  First,  the  eight  species 
production  equations  are  solved  simultaneously,  followed  by  the 
separate  vibrational  relaxation  equations,  and  then  the  electron 
energy  equation.  After  these  equations  are  solved  for  the  increment 
h,  the  integrated  equations  for  conservation  of  mass,  momentum,  and 
energy  along  with  the  equation  of  state  must  be  satisfied.  This  is 
accomplished  by  iteration.  Newton's  method  is  used. 

A stability  and  truncation  error  analysis  has  not  been  performed 
for  the  numerical  integration  method  used  here;  however,  SHFRMAN  (1968) 
does  prove  convergence  for  use  of  the  simple  approximate  form.  The 
corresponding  computer  program  is  not  in  the  production  stage;  all 
input  data  do  not  produce  stable  solutions  and  sufficient  cases  have 
not  been  calculated  to  establish  more  than  just  general  guidelines 
for  choices  of  initial  conditions . 


All  calculations  were  performed  on  a CDC  6500  digital  computer 
and  the  program  is  coded  in  FORTRAN  IV.  Computer  program  documenta- 
tion, brief  subroutine  descriptions,  and  input  data  information 
(initial  conditions)  is  not  provided  in  this  report  but  can  be  made 
available  upon  request  to  the  author. 
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4 . RESULTS 

4 . 1 Post-Shock  Equilibrium  Conditions 


The  shocked  equilibrium  conditions  attained  following 
relaxation  using  the  airshock  relaxation  code  have  been  compared  with 
calculations  from  WILLETT  and  LEHTO  (1958)  for  a Mach  number  range 
M = 7-36  with  ambient  conditions:  T = 300°K,  p = 5 Torr.  The  results 

of  WILLETT  and  LEHTO  (1958)  are  obtained  by  solving  the  normal  shock 
(Rankine-Hugoniot)  relations  without  considering  relaxation  mechanisms. 
The  corresponding  properties  of  the  final  equilibrium  state  agree 
within  1%  for  the  variables  of  temperature,  density,  pressure,  and 
velocity.  Electron  number  density  results  obtained  with  the  airshock 
relaxation  code  agree  on  the  order  of  15%  with  equilibrium  properties 
tabulated  by  HILSENRATH  and  KLEIN  (1965) . Although  the  detailed 
results  are  not  presented  here,  the  comparisons  do  indicate  that  proper 
post-shock  equilibrium  conditions  are  calculated  with  the  airshock 
relaxation  code. 

4 . 2 Relaxation  Region  Solution 

Figures  1-4  display  the  relaxation  region  solution  obtained 
for  the  flow  conditions:  M = 16.47,  p = 5 Torr,  T = 300°K,  which  are 

in  the  regime  studied  by  LIN  and  TEARE  (1963)  and  SCHAFER  and  FROHN 
(1972) . LIN  and  TEARE  (1963)  employ  "one-way  coupling"  between  the 
processes  of  dissociation  and  ionization  which  is  satisfactory  when 
the  degree  of  ionization  is  less  than  'vli;  the  temperature,  density, 
and  chemical  composition  profiles  are  determined  (including  O2  and  N2 
vibrational  relaxation)  without  considering  ionization  processes.  The 
ionization  solution  is  then  completed  using  the  determined  flow  field. 
SCHAFER  and  FROHN  (1972)  solve  the  coupled  chemical  and  ionizational 
relaxation  equations  without  considering  vibrational  relaxation. 

A comparison  between  Figure  1 and  Figure  3 indicates  that 
N0+  is  the  predominant  ion.  Figure  2 displays  the  early  disappearance 
of  02  and  the  appearance  of  an  NO  peak.  These  are  expected  results 
for  these  flow  conditions  (LIN  and  TEARE  (1963)) . The  temperature 
profiles  presented  in  Figure  4 show  strong  coupling  between  N2  vibra- 
tional energy  and  the  electron  energy.  This  is  in  agreement  with  the 
f ree-electron/N2-vibrational  model  and  results  discussed  by  ZHELEZNYAK 
and  MNATSAKANYAN  (1968) . Initially,  for  the  present  calculations,  the 
free  electrons  are  assumed  to  be  in  equilibrium  with  the  N2  vibrators. 
At  the  point  where  the  particle  time  increment  between  calculation 
steps  becomes  an  order  of  magnitude  larger  than  the  time  required 
for  free  electrons  to  establish  a Maxwell -Boltzmann  distribution, 
separate  ^-vibrator  and  free-electron/bound-electron  energy  relaxa- 
tion is  permitted.  This  transition  is  indicated  in  Figure  4 by  the 
appearance  of  a dip  in  the  electron  temperature.  After  several  steps 
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(10  for  this  calculation) , the  electrons  approach  equilibrium  with 
the  N2-vibrators . 

The  o2  vibrational  temperature  profile  in  Figure  4 exhibits 
a pronounced  peak  whereas  the  No  vibrational  temperature  does  not  for 
these  flow  conditions.  This  behavior  is  described  in  the  paper  by 
TREANOR  and  MARRONE  (1962)  from  which  the  model  for  vibrational 
relaxation  (with  no  N2-vibrational/free-electron  coupling)  used  in 
these  calculations  is  taken. 


Figures  5 and  6 correspond  to  the  relaxation  region  solution 
for  the  conditions:  M = 26.76,  p = 5 Torr,  T = 300°K.  The  electron 

number  density  profile  in  Figure  5 is  quite  different  from  the  profile 
in  Figure  1 (M  = 16.47).  Figure  6 indicates  that  N2+  and  0+  are  the 
predominant  ions.  The  associative  ionization  reactions,  N+o£NO+  + e 
is  the  main  ionization  mechanism  up  to  Mv26 ; however,  for  M>26 
reactions  such  as  N+N?N^+  + e and  0+e20+  + 2e  predominate.  The  general 
trend  is  that  the  main  ionization  mechanism  shifts  from  associative 
ionization  to  electron-neutral  collisional  ionization  for  increasing 
Mach  numbers  beyond  Mv26. 


BIBERMAN  and  YAKUBOV  (1965)  and  ZHELEZNYAK  and  MNATSAKANYAN 
(1968)  provide  additional  discussion  on  ionization  mechanisms  for 
M>26 . * The  former  present  an  ionizational  relaxation  model  appropri- 
ate for  M>29;  the  lower  limit  due  to  the  restrictive  assumption  that 
the  dissociation  process  for  02  and  N2  are  considered  complete  before 
the  onset  of  ionization.  ZHELEZNYAK  and  MNATSAKANYAN  (1968)  provide 
an  ionizational  relaxation  model  satisfactory  for  M>26.  In  this  case 
the  dissociation  of  N2  is  solved  simultaneously  with  ionization  and 
C>2  is  assumed  completely  dissociated,  initially. 

4.3  Ionizational  Relaxation 


Figure  7 provides  a comparison  between  available  experimental 
data  and  the  present  calculations  for  ionizational  relaxation  time  as 
a function  of  Mach  number  for  the  range  M = 7-36.  Other  solutions  for 
portions  of  this  Mach  number  range  which  have  been  discussed  briefly 
in  section  4.2  Relaxation  Region  Solution  of  text  are  presented  in 
BIBERMAN,  MNATSAKANYAN,  and  YAKUBOV  (1971).  Experimental  data  from 
eight  sources  are  included  in  Figure  7.**  Ambient  pressures  for  these 
experiments  range  from  0.05  to  10  Torr.  Since  the  ionization  mech- 
anism is  a binary  collision  process,  ionizational  relaxation  times 
are  pressure-scaled  for  purposes  of  data  reduction.  Pressure  scaling 
has  been  verified  for  these  calculations  which  do  not  include  effects 
of  radiation. 


♦ZHELEZNYAK  (1973)  discusses  ionization  relaxation  in  nitrogen  for 
M = 50-73. 

♦♦Experimental  data  plotted  in  Figure  7 are  taken  from  graphs  provided 
by  the  cited  references. 
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The  Mach  number  range  presented  in  Figure  7 can  be  separated 
into  three  subregions  established  by  the  distribution  of  exoerimental 
data.  The  first  subregion,  M = 7-20,  is  well  defined  by  data  of 
SCHAFER  and  FROHN  (1972)  and  LIN,  NEAL,  and  FYFE  (1962).  Other  data 
found  in  the  literature  are  included  in  the  figure.*  For  this  region 
the  definition  of  ionizational  relaxation  time  (laboratory  coordinates) 
for  the  theoretical  model  is  taken  to  be  the  time  for  the  electron 
number  density  to  attain  90%  of  the  peak  value.  The  calculated  curve 
compares  best  with  the  data  of  SCHAFER  and  FROHN  (1972)  for  the  range 
M = 7-20. 


No  data  have  been  located  in  the  range  M = 20-26  for 
comparisons  with  calculations.  Two  additional  definitions  for  ioniza- 
tional relaxation  time  are  employed  for  the  computations  which  give 
equivalent  results  for  this  region  with  agreement  on  the  order  of  10%. 
The  definition  used  in  Figure  7 is  as  follows:  The  time  (laboratory 

coordinates)  between  the  shock  front  and  the  time  at  which  a line 
drawn  through  the  maximum  slope  of  the  electron  number  density  profile 
reaches  the  equilibrium  value.  Another  definition  which  is  also 
checked  is  one  used  by  WILSON  (1966)  for  analysis  of  experimental 
results  in  the  range  M = 26-36.  An  infrared  radiation  model  which 
includes  bremsstrahlung  and  recombination  continuum  for  oxygen  and 
nitrogen  (from  Wilson  (1966))  is  used  to  calculate  a radiation  profile 
for  the  relaxation  region  at  a wavelength  of  6.1y.  The  relaxation 
time  is  then  defined  to  be  the  time  (laboratory  coordinates)  between 
the  shock  front  and  the  time  at  which  a line  drawn  through  the 
maximum  slope  of  the  infrared  radiation  profile  reaches  the  equilibrium 
value . 


The  range  M = 20-26  contains  a portion  of  the  transition 
region  for  the  shift  in  the  ionizational  relaxation  mechanisms  dis- 
cussed in  section  4.2  Relaxation  Region  Solution  of  text.  The  cal- 
culations evidence  this  shift  in  Figure  7.  The  lack  of  data  in  this 
region  does  leave  to  question  the  proper  location  for  the  shift  in 
ionization  mechanisms. 

The  comparison  between  the  present  calculations  and  the 
experimental  results  of  WILSON  (1966)  for  the  range  M = 26-36  indi- 
cates only  fair  qualitative  agreement.  The  definition  of  ionizational 
relaxation  time  for  this  region  is  the  same  as  that  for  M = 20-26  and 
is  equivalent  to  the  definition  used  by  WILSON  (1966) . Because  of  the 
transition  in  ionization  mechanisims  there  is  a peak  in  the  relaxation 
times  as  a function  of  Mach  number.  The  theoretical  and  experimental 
locations  for  this  peak  do  not  agree.  Modifications  have  not  been 
made  to  the  model  which  account  for  this  disagreement.  No  radiation 
model  has  been  included  for  iteration  with  the  flow  field  model  for 
strong  radiative-gasdynamic  coupling  effects  described  by  CARLSON 
(1972)  for  M>29.  It  should  be  pointed  out,  however,  that  ZHELEZNYAK 
and  MNATSAKANYAN  (1968)  obtain  good  agreement  with  the  data  of 

*LIN,  NEAL,  and  FYFE  (1962)  discuss  the  discrepancies  between  their 
data  and  the  earlier  results  of  NIBLETT  and  BLACKMAN  (1958)  and  the 
data  of  MANHEIMER-TIMNAT  and  LOW  (1959). 
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WILSON  (1966)  for  these  flow  conditions  without  including  effects  of 
radiative-gasdynamic  coupling.  Also,  no  adjustments  of  selected 
reaction  rates  nor  addition  of  kinetic  equations  for  individual 
excited  states  of  atomic  nitrogen  and  oxygen  (as  considered  by 
ZHELEZNYAK  and  MNATSAKANYAN  (1968))  have  been  included  to  obtain 
better  agreement  between  the  results  of  the  present  formulation  and 
the  data  of  WILSON  (1966) . 

4.4  O2  Dissociation 


The  variation  of  O2  half-dissociation  length  with  Mach  number 
for  the  range  M = 8-22  is  given  in  Figure  8.  The  solid  curve  shown  in 
the  figure  is  derived  from  results  obtained  with  the  airshock  relaxa- 
tion computer  code.  For  these  calculations  the  O2  half-dissociation 
length  is  defined  as  the  location  in  the  relaxation  region  where  the 
mole  fraction  of  O2  attains  the  value  Xo2= ( (Xc^) INITIAL  + (XC^) FINAL^ • 

Experimental  results  from  STRICKER  and  LOW  (1972)  and  LIN  and  FYFE 
(1961)  are  included  for  comparison.  Pressure  scaling  is  used  to 
convert  the  experimental  data  to  the  proper  ambient  conditions.  The 
calculations  are  on  the  high  side  but  are  within  a factor  of  two  for 
the  low  Mach  number  region  M = 8-10  and  the  calculations  are  within 
the  scatter  of  data  for  the  high  Mach  number  region,  M = 14-21.  The 
calculations  for  the  latter  range,  however,  do  tend  to  fall  off  more 
rapidly  than  the  data  indicate.  STRICKER  and  LOW  (1972)  present 
additional  theoretical  curves  for  O2  half-dissociation  length 
dependence  on  Mach  number  which  do  not  agree  as  well  with  the 
experimental  results  for  the  entire  range  as  do  the  present 
calculations . 

4.5  Vibrational  Relaxation 


Results  for  vibrational  relaxation  of  and  N2  using  the 
airshock  relaxation  computer  code  are  given  in  Figure  9 . Vibrational 
relaxation  in  particle  coordinates,  not  laboratory  coordinates  as 
used  in  Figures  1-7,  is  plotted  versus  temperature. 

The  expression  for  vibrational  relaxation  is  defined  in 
APPENDIX  B.  The  main  contribution  to  vibrational  relaxation  and  in 
many  cases  the  only  contribution  considered,  is  collisional-vibrational 
energy  exchange.  This  contribution  pressure  scales  due  to  collisional 
dependence  thus  justifying  the  use  of  pi  ordinates  in  Figure  9 . The 
Landau-Teller  theory  (MILLIKAN  and  WHITE  (1963))  postulates  inverse 
cube  root  temperature  dependence  for  vibrational  relaxation  which  is 
the  reason  behind  the  choice  of  the  temperature  function  for  the 
abscissa  of  Figure  9.  Also  included  in  the  figure  are  empirical  fits 
to  experimental  data  for  vibrational  relaxation  of  O2  and  N2  estab- 
lished by  MILLIKAN  and  WHITE  (196  3)  usincr  the  Landau-Teller  theory  as 
a guide  for  functional  dependence  of  important  parameters  involved 
with  the  interaction.  The  MILLIKAN  and  WHITE  (1963)  model  only 
considers  collisional-vibrational  energy  exchange.  This  same  model 
is  used  for  the  collisional-vibrational  energy  exchanae  contribution 
in  the  present  calculation.  For  N2 , additional  contributions  to 
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vibrational  relaxation  from  the  dissociation/association  process  for 
N2  and  free-electron/N2~vibrator  energy  exchange  are  included.  For 
O2,  the  additional  contribution  from  the  dissociation/association 
process  for  O2  is  included.  The  difference  between  the  MILLIKAN  and 
WHITE  (1963)  results  and  the  present  results  in  Figure  9 is  attrib- 
utable in  most  part  to  the  additional  contributions  discussed  above 
for  the  present  model.  Another  reason  for  differences  between  the 
two  results  follows  from  the  definitions  of  pressure  and  temperature 

-1 

used  to  convert  the  present  calculations  to  px-T  coordinates*.  The 
local  temperature  and  pressure  values  used  for  the  airshock  relaxation 
results  correspond  to  the  flow  conditions  just  downstream  of  the 
translational/rotational  shock;  that  is,  the  initial  temperature  and 
pressure  in  the  relaxation  region  before  chemical,  vibrational, 
electronic,  and  ionizational  relaxation  processes  have  begun.  However, 
if  the  airshock  relaxation  results  in  Figure  9 were  corrected  to 
correspond  to  the  local  pressure  and  temperature  at  the  location  of 
the  peak  vibrational  temperature , which  would  represent  the  end  of  the 
vibrational  relaxation  region,  then  the  results  would  be  modified  at 
most  on  the  order  of  10%;  approximately  7%  change  for  temperature  and 
3%  change  for  pressure.  It  is  for  this  reason  that  the  results  are 

-1 

. 3 

plotted  using  px-T  coordinates  along  with  the  MILLIKAN  and  WHITE 
(1963)  results  even  though  contributions  from  relaxation  processes 
other  than  direct  collisional-vibrational  energy  exchange  are 
included. 


The  net  effect  of  including  the  contributions  of  the 
dissociation/association  processes,  and  in  the  case  of  N2  also 
including  free-electron/N2~vibrator  coupling,  in  the  vibrational 
relaxation  model  is  to  reduce  the  relaxation  time  as  evidenced  in 
Figure  9.  The  higher  the  Mach  number,  or  temperature,  the  greater 
the  reduction  in  relaxation  time  for  both  N2  and  O2 . 

For  the  higher  Mach  numbers,  M>26  or  (T-^/^)  values  less 
than  0.030  in  Figure  9,  the  vibrational  relaxation  model  used  is  not 
satisfactory.  The  px  range  is  of  the  order  of  collision  times  and 
assumptions  for  the  vibrational  relaxation  model  such  as  the  estab- 
lishment of  Maxwell -Boltzmann  distributions  are  just  not  satisfied. 
However,  in  this  Mach  number  range,  M>26,  vibrational  contributions 
to  other  relaxation  processes  such  as  ionization  become  less 
important.  The  calculations  indicate  that  the  dissociation  of  O2  can 
be  considered  complete  before  the  onset  of  ionization  for  M>26 
(ZHELEZNYAK  and  MNATSAKANYAN  (1968))  and  N2  can  be  considered  comple- 
tely dissociated  before  ionization  begins  for  M>29  (BIBERMAN  and 
YAKUBOV  (1965)).  Therefore,  only  improvement  of  the  vibrational 
relaxation  model  for  N2  through  the  range  M = 26-29  is  needed. 


* ^1 

3 

px  versus  T 
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4 . 6 Electron  Energy  Relaxation 

An  example  of  an  electron  temperature  profile  for  the 
relaxation  region  is  provided  in  Figure  4.  The  discussion  presented 
in  section  4.2  Relaxation  Region  Solution  of  text  will  not  be  repeated 
here  except  to  again  point  out  the  strong  energy  coupling  between  the 
N2  vibrators  and  the  free  electrons  indicated  in  the  figure.  The 
electron  energy  equation  which  is  used  to  calculate  this  profile  is 
defined  in  APPENDIX  C. 

For  the  lower  Mach  number  range  M = 7-12*  calculations 
indicate  that  the  electron  energy  approaches  final  shocked  equilibrium 
conditions  very  rapidly,  on  the  order  of  the  time  it  takes  to  establish 
a separate  electron  Maxwell-Boltzmann  distribution.  The  electron 
temperature  then  remains  essentially  constant  throughout  the  rest  of 
the  relaxation  region.  Therefore,  calculation  time  can  be  reduced  by 
assuming  that  the  electron  temperature  equals  the  post-shock  equili- 
brium temperature  throughout  the  relaxation  region  for  the  above 
stated  flow  conditions. 

Within  the  range  M = 12-20  the  main  contribution  to  the 
calculation  by  including  a separate  electron  temperature  is  in  reduc- 
ing the  amount  of  energy  transferred  to  the  N2  vibrators  over  most  of 
the  relaxation  region.  When  Te  is  constrained  to  follow  the  atom 
translational  temperature,  then  the  N2  vibrators  appear  to  be  in 
near  equilibrium  with  translation  for  most  of  the  relaxation  region. 
However,  if  a separate  electron  temperature  is  permitted,  then  the 
free  electrons  and  the  N2  vibrators  are  in  near  equilibrium  for  a 
major  portion  of  the  relaxation  region  and  they  together  approach 
equilibrium  with  translation  much  more  slowly.** 

Above  M^20  the  main  contribution  for  including  a separate 
electron  temperature  shifts  from  modifying  the  ^-vibrator  energy 
distribution  to  proper  evaluation  of  reaction  rates  for  electron- 
atom  collisional  ionization  and  electron-ion  recombination.  These 
reaction  rates  depend  on  the  free  electron  energy  distribution  (tem- 
perature) and  an  evaluation  of  the  reaction  rates  at  the  higher  atom 
translational  temperature  is  in  error. 

For  M>25  considerable  calculation  time  is  consumed  by 
considering  a separate  electron  temperature.  This  calculation  time 
can  be  greatly  reduced  by  setting  the  electron  temperature  equal  to 
the  N2  vibrational  temperature  at  all  points  in  the  relaxation  region. 
This  is  a reasonable  assumption  because  of  the  strong  energy  coupling 
between  the  free  electrons  and  the  N2  vibrators?  there  is  near 
equilibrium  between  them  over  most  of  the  relaxation  region. 


*Mach  number  ranges  quoted  in  this  section  are  only  approximate.  A 
sufficient  number  of  cases  have  not  been  investigated  to  make  truly 
accurate  determinations . 

**The  free-electron/N2~vibrator  coupling  terms  are  not  considered 
when  the  free  electron  temperature  is  constrained  to  follow  the 
translational  temperature. 
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5.  SUMMARY  AND  RECOMMENDATIONS 

A model  for  the  relaxation  region  behind  shock  waves  in  air  has 
been  developed.  It  is  not  the  first  such  model.  However,  other 
models  to  date  are  only  satisfactory  for  portions  of  the  Mach  number 
range,  M = 7-36,  considered  here.  Within  this  range  the  ionizational 
mechanism  shifts  from  associative  ionization  to  electron-neutral 
collisional  ionization  and  many  air  chemical  reactions  change  in 
relative  importance  with  respect  to  one  another. 

This  model  has  been  developed  for  use  in  determining  radiation 
output  from  airshocks . The  airshock  relaxation  region  is  determined 
so  that  the  flow  properties  may  be  used  to  calculate  the  radiation 
properties  for  the  region.  One-way  coupling  is  employed.  The  flow 
field  determines  the  radiation  field  but  it  is  assumed  that  the 
radiation  field  does  not  affect  the  flow  field  appreciably.  The 
present  calculations  may  be  considered  the  first  step  in  an  iteration 
process  in  which  the  computations  are  terminated  when  the  flow  field 
and  radiation  field  results  for  consecutive  cycles  agree  within  stated 
limits.  The  radiation  calculations  have  not  been  performed  using  this 
model . 

The  summary  of  results  to  follow  in  subsequent  paragraphs  is 
separated  into  five  Mach  number  ranges  which  have  been  determined  by: 

1.  The  distribution  of  experimental  data. 

2.  The  shift  of  the  predominant  ionization  mechanism  from 
associative  ionization  to  electron-neutral  collisional  ionization. 

3.  The  change  in  character  of  the  separate  electron/atom  energy 
distributions . 

In  addition  to  comments  concerning  agreement  between  calculated  and 
experimental  results  for  each  Mach  number  range,  important  considera- 
tions such  as  the  predominant  ionization  mechanism,  the  status  of  O2 
and  N2  dissociation  and  vibrational  relaxation,  and  the  properties  of 
the  separate  electron/atom  energy  distributions  are  noted.  Refer 
to  Figures  7,  8,  and  9 for  comparisons  concerned  with  ionizational 
relaxation  times,  C>2  dissociation  distances,  and  vibrational 
relaxation  times  for  O2  and  N2  respectively. 

M = 7-12  Associative  ionization  is  the  predominant  ionization 
mechanism. 

C>2  and  N2  dissociation  and  vibrational  relaxation  must  be 
considered. 

Electron  temperature  is  essentially  equal  to  the  post-shock 
equilibrium  temperature  throughout  the  relaxation  region. 
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Calculations  for  post-shock  equilibrium  conditions  agree 
with  standard  Rankine-Hugoniot  results. 

Calculations  for  ionizational  relaxation  times  are  within 
scatter  of  data. 

Calculations  for  O2  dissociation  distances  are  within  scatter 
of  data. 

Reasonable  values  for  O2  and  N2  vibrational  relaxation  times 
are  calculated  with  model. 

M = 12-20  Associative  ionization  is  the  predominant  ionization 
mechanism. 

O2  and  N2  dissociation  and  vibrational  relaxation  must  be 
considered. 

Separate  electron  and  N2~vibrator  energy  distributions  are 
in  near  equilibrium  for  most  of  the  relaxation  region. 

Calculations  for  post-shock  equilibrium  conditions  agree 
with  standard  Rankine-Hugoniot  results. 

Calculations  for  ionizational  relaxation  times  are  within 
scatter  of  data. 

Calculations  for  O2  dissociation  distances  are  within  scatter 
of  data.  No  data  available  above  M = 16. 

Reasonable  values  for  O2  and  N2  vibrational  relaxation  times 
are  calculated  with  model. 

M = 20-26  Associative  ionization  is  the  predominant  ionization 
mechanism. 

C>2  and  N2  dissociation  and  vibrational  relaxation  must  be 
considered . 

Separate  electron  and  N2~vibrator  energy  distributions 

are  in  near  equilibrium  throughout  most  of  the  relaxation 
region . 

Calculations  for  post-shock  equilibrium  conditions  agree 
with  standard  Rankine-Hugoniot  results. 

No  experimental  data  for  ionizational  relaxation  times  are 
available.  Theoretical  model  indicates  abrupt  changes 
in  Mach  number  dependance  of  ionizational  relaxation  times 
within  range  M = 25-26. 

Reasonable  values  for  O2  and  N2  vibrational  relaxation  times 
are  calculated  with  model. 

M = 26-29  The  predominant  ionization  mechanism  shifts  from  associative 
ionization  to  electron-neutral  collisional  ionization. 

N2  dissociation  and  vibrational  relaxation  must  be 
considered. 

Electron  temperature  is  essentially  equal  to  local  N2 

vibrational  temperature  throughout  most  of  the  relaxation 
region . 

Calculations  for  post-shock  equilibrium  conditions  agree 
with  standard  Rankine-Hugoniot  results. 
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Calculations  for  ionizational  relaxation  times  follow 

qualitatively  the  abrupt  chanqe  in  Mach  number  dependence 
of  data.  However,  the  calculations  do  fall  below  the  data 
and  the  pronounced  peak  indicated  in  the  data  (Figure  7) 
is  not  followed  by  the  calculations. 

The  vibrational  relaxation  model  is  not  adequate  for  the 
range.  Unrealistically  short  relaxation  times  are 
calculated  for  N2  vibrators . 

M = 29-36  Electron-neutral  collisional  ionization  is  the  predominant 
ionization  mechanism. 

N2  vibrational  relaxation  must  be  considered. 

Electron  temperature  is  essentially  equal  to  local  N2 

vibrational  temperature  throughout  most  of  the  relaxation 
region . 

Calculations  for  post-shock  equilibrium  conditions  agree 
with  standard  Rankine-Hugoniot  results. 

Calculations  for  ionizational  relaxation  times  do  not  follow 
Mach  number  dependence  of  data. 

The  vibrational  relaxation  model  is  not  adequate  for  this 
range.  Unrealistically  short  relaxation  times  are 
calculated  for  N2  vibrators. 

The  following  modifications  to  the  airshock  relaxation  model  are 
recommended  so  as  to  obtain  better  agreement  between  calculations 
and  experimental  data  for  the  Mach  range  M = 26-36. 

1.  The  existing  set  of  chemical  reaction  rates  should  be  adjusted 
on  an  individual  reaction  basis  in  order  to  optimize  aareement  between 
calculations  and  data. 

2.  Reactions  which  account  for  individual  excited  state 
populations  of  atomic  nitrogen  and  oxyoen  should  be  added  to  the  set 
of  reactions.* 

The  above  two  modifications  to  the  model  should  then  be  followed 
by  iterations  between  flow  field  calculations  (the  type  of  calculations 
presented  in  this  report)  and  radiation  output  calculations  (to  be 
developed) . The  actual  effects  of  radiation  of  the  flow  field  can 
then  be  determined. 

Also,  the  vibrational  relaxation  model  needs  to  be  modified  for 
the  same  Mach  number  range,  M = 26-36.  The  present  model  defines 
vibrational  temperatures  in  the  relaxation  region  at  unrealistically 
early  times.  This  procedure  can  be  corrected  by  not  defining  a 


*ZHELEZNYAK  and  MNATSAKANYAN  (1968)  obtained  good  agreement  with 
experimental  data  for  ionizational  relaxation  times  within  this 
Mach  number  range  by  having  included  these  reactions  and  by  not 
having  included  radiative-gasdynamic  coupling. 
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vibrational  temperature  before  a sufficient  number  of  collisions  have 
occurred  to  establish  a Maxwell-Boltzmann  distribution  of  vibrators. 
This  modification  will  also  affect  the  early-time  definition  of  the 
electron  temperature.  The  vibrational  temperatures  (N2 ) and  the 
electron  temperature  can  be  evaluated  at  the  pre-shock  (ambient) 
conditions  throughout  this  initial  relaxation  region. 

NOTE  ADDED  IN  PROOF 

Two  recent  papers  concerned  with  results  of  interest  to  the 
present  investigation  came  to  the  attention  of  the  author 
prior  to  publication  of  this  report. 

Chan  and  Glass*  present  experimental  data  which  fall  within 
the  data  boundaries  of  WILSON  (1966)  and  follow  closely  the 
ionizational  relaxation  theory  of  ZHELEZNYAK  and  MNATSAKANYAN 
(1968)  for  the  Mach  range  M = 28-35.  These  results  provide 
an  additional  indication  that  the  airshock  relaxation  model, 
the  subject  of  this  report,  should  be  modified  as  discussed 
in  section  5.  SUMMARY  AND  RECOMMENDATIONS  of  the  text  to 
obtain  better  agreement  with  data  in  this  flow  regime. 

Lazdinis  and  Petrie**  report  on  recent  theoretical  and 
experimental  results  concerned  with  free  electron  and 
^-vibration  energy  coupling  for  weakly  ionized  nozzle 
expansions  of  shock-heated  nitrogen.  Although  the  nozzle 
expansion  flow  conditions  are  quite  different  from  airshock 
relaxation  flow  conditions , the  formulations  for  the  .electron 
energy  equations  are  similar;  differences  arise  due  to  the 
relative  importance  of  the  individual  electron  energy  flux 
terms.  The  main  coupling  mechanism  for  both  flows  is  the 
same;  that  is  free-electron/^-vibrator  energy  exchange. 


*Chan,  S.  K.  and  Glass,  I.  I.,  "Radiative  Relaxation  Behind  High-Speed 
Shock  Waves  in  Air,"  Physics  of  Fluids  1T_,  No.  4,  19  74 

**Lazdinis,  S.  S.  and  Petrie,  S.  L. , "Free  Electron  and  Vibrational 
Temperature  Nonequilibrium  in  High  Temperature  Nitrogen" , Physics 
of  Fluids  17,  No.  8,  1974 
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6.  LIST  OF  SYMBOLS 

A 

- Constant  used  in  definition  of  chemical  reaction 
rate,  units  are  (CM3)D/SEC  for  which  j = 0,1,  or  2 
depends  on  the  reaction  being  of  first,  second,  or 
third  order. 

H 

- Coefficient  of  constant  term  in  the  approximate  form 
for  f^. 

Aij 

- Coefficient  matrix  used  for  system  of  species 
conservation  equations. 

ae 

- Average  electronic  excitation  cross  section  for 
free-electron  impact,  CM^. 

B 

- Constant  used  in  definition  of  chemical  reaction  rate 
for  temperature  power  dependence. 

C 

- Constant  used  in  definition  of  chemical  reaction  rate 
which  represents  activation  energy  in  temperature 
units,  °K . 

i 

- Net  production  rate,  forward  minus  reverse,  for  the 
it*}  electron-neutral  collisional  ionization  reaction 
CM3/SEC . 

Ei 

- Electron  excitation  energy  for  chemical  species  i, 
ERGS . 

Ev 

- Vibrational  energy,  ERGS. 

(Ev) VIB. 

E0  - Vibrational  energy  evaluated  at  local  translational/ 

vibrational  equilibrium  conditions,  ERGS. 

EIi 
e (t) 

- Ionization  energy  of  chemical  species  i,  ERGS. 

- Partition  function  for  N2  harmonic  oscillator. 

fi 

- Derivative  of  species  y^  with  respect  to  x. 

fij 

- Derivative  of  species  y^  with  respect  to  x evaluated 
at  the  jth  point. 

9i 

- Statistical  weight  of  i*-*1  electronic  energy  level . 

h 

- Step  size  for  numerical  integration,  CM. 
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h(T) 

- Enthalpy,  ERGS/GM 

hF 

- Enthalpy  of  formation,  ERGS/GM-MOLE . 

K 

- Chemical  equilibrium  constant 

Ki 

- Rate  of  production  of  species  i per  unit  volume, 
CM“3/SEC. 

k 

- Boltzmann  constant,  ERG/°K. 

kd 

- Dissociation  reaction  rate  corrected  for  vibrational 
relaxation  effects,  CM3/SEC. 

(kd)viB.  EQ. 

- Dissociation  reaction  rate  evaluated  for  translational/ 
vibrational  equilibrium,  CM3/SEC. 

kf,  kr 

- Forward  and  reverse  chemical  reaction  rates,  units 
are  defined  for  SYMBOL  A. 

(kAI^  f ' ^kAI ^ r ‘ 

- Forward  and  reverse  associative  ionization  reaction 
rates , CM3/SEC . 

kRR 

- Chemical  reaction  rate,  units  are  defined  for 
SYMBOL  A. 

M 

- Mach  number. 

mi 

- Mass  of  ion  species  i,  GM. 

N 

- Number  of  discrete  energy  levels  below  dissociation 
energy  for  cut-off  harmonic  oscillator  models  of 
O2  and  N2. 

n 

. . — ? 

- Total  particle  number  density,  CM  . 

ne 

- Electron  number  density,  CM-3. 

Pij 

* Coefficient  matrix  in  f^  approximate  form  evaluated 
for  fi  and  y j . 

pe 

- Average  rate  constant  for  electronic  excitation  by 
free-electron  impact,  CM3/SEC. 

P 

- Static  pressure,  DYNE/CM2  (also  Torr) . 

Qi 

- Electron  energy  flux  terms  defined  in  Appendix  C, 
i=l , 5 ; GM/SEC3-CM. 

R 

- Mass  gas  constant,  ERG/°K-GM. 

T 

- Temperature,  translational  temperature,  °K. 
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Tv,TVN2,TV02 

U 

u 


Xi 

x 

y 

y± 

Yij 
Yeq 
Zl,  Z2 


'v 


P 

T 


tl 

TP 

SUBSCRIPTS 


i, j ,k,  l,m 


o 


V 


- Electron  temperature,  °K. 

- (Tv-1  + T"1)"1  °K . 

- Vibrational  temperature,  N2  and  O2 , °K. 

- Velocity  of  shock  front  into  undisturbed  gas,  CM/SEC. 

- Gas  velocity  behind  shock  wave  in  the  coordinate 
system  fixed  to  wave  front,  CM/SEC. 

- Electron  mean  velocity,  CM/SEC . 

- Mole  fraction  of  species  i. 

- Distance  coordinate,  independent  variable,  CM. 

- Dependent  variable,  function  of  x. 

- Dependent  variable  for  soecies  i,  function  of  x. 

- Dependent  variable  for  species  i evaluated  at  point 
j,  function  of  x. 

- Dependent  variable  evaluated  at  local  equilibrium 
conditions,  function  of  x. 

- Uncertainty  parameters  used  in  chemical  reaction 
rate  definition. 

- Characteristic  vibrational  temperature,  °K. 

- Characteristic  electronic  temperature  for  energy 
level  i,  °K. 

- Gas  density,  GM/CM^ . 

- Relaxation  time,  ySEC. 

- Relaxation  time,  laboratory  coordinates,  ySEC. 

- Relaxation  time,  particle  coordinates,  ySEC. 

- Dummy  subscripts . 

- Free  stream  conditions,  pQ,  pQ,  TQ,  (Xi)Q. 

- Denotes  vibrational  relaxation  quantity. 
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APPENDIX  A:  CHEMICAL  REACTION  RATES 

The  chemical  reaction  rates  used  in  the  airshock  chemistry  model 
are  listed  at  the  end  of  this  appendix  in  Table  A-l . Types  of  reac- 
tions include  dissociation,  rearrangement,  ionization,  and  positive- 
ion  charge  exchange  with  their  respective  inverse  reactions.  The 
reaction  rates  are  taken  from  three  sources:  reference  A1  - BORTNER 

(1963),  reference  A2  - BORTNER  (1969),  and  reference  A3  - BORTNER  and 
BAURER  (1972).  The  rates  are  listed  in  either  the  forward  (kf)  or 
reverse  (kr)  direction  and  are  presented  in  the  form* 

kRR  = (A+Zl)  (T/300)b±z2  exp  (-C/T)  , 

with  units  (CM3)3  SEC-'1'  for  which  j = 0,  1,  or  2 depends  on  the 
reaction  being  of  first,  second,  or  third  order,  respectively.  Z1 
and  Z2  represent  the  uncertainty  estimates.  Better  agreement  between 
the  present  calculations  and  cited  experimental  results  for  ioniza- 
tional  relaxation  times  is  obtained  using  the  rate  parameters 
A'  = A-Zl  and  B'  = B-Z2. 

The  equilibrium  constant  is  used  to  relate  the  rates  of  the 
opposing  directions  by 


K 


Actually  this  relation  only  holds  for  reactions  in  which  the  same 
states  (electronic,  vibrational,  and  rotational)  are  attained  for  pro- 
ducts and  reactants.  Considering  the  available  data  and  the  corre- 
sponding uncertainty  limits , the  use  of  the  above  relation  is 
reasonable.  However,  for  the  dissociation  reactions  of  O2  and  N2  the 
vibration-dissociation  coupling  factor  proposed  by  HAHMERLING , TEARE 
and  KIVEL  (1959)  is  used  to  correct  the  reaction  rates  for  the 
dissociation  process . 


k.  = 

d d VIB.EQ. 


-NSv/Tm  Sv/T„ 
(1-e  m)  (e  -1) 


. ®v/Tm  . (e®v/T-l) 

N(e  -1) 


The  dissociation  reaction  rates  are  modified  to  account  for  the 
distribution  of  the  vibrational  states  being  characteristic  of  the 
local  vibrational  temperature  and  not  the  local  translational 


*See  6.  LIST  OF  SYMBOLS  of  text. 
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temperature.  The  assumption  that  Tv  = T when  actually  Tv  ^ T 
incorrectly  prescribes  the  distribution  of  the  vibrational ly,  excited 
states.  This  affects  the  dissociation  reaction  rate  since  the  highly 
excited  vibrational  states  are  more  easily  dissociated. 

The  ionization  reaction  rate  dependence  on  electron  temperature 
is  outlined  in  2.  THEORY  of  text. 

Fifty  reactions  are  considered,  though  only  eight  independent 
equilibrium  constants  are  required.  The  eight  reactions  selected  are 

0 2 t 0 + 0 

N2  t N + N 

NO  t N + 0 

N t N+  + e 

0 t 0+  + e 

NO  t N0+  + e 
02  t 0£  + e 

N2  t + e 

and  the  corresponding  species  conservation  equations  may  be  expressed 
in  the  form  (SHERMAN  (1968)) 


Xq 

CM 

1 

o 

X 

Xq 

a 

3xo 

3xo 

2X0 

2X0 

o 

X 

CM 

Xq 

*n 

*N 

(XN 

0 

3xn 

3xn 

2xn 

2xn 

2xn 

Xn 

XN0 

XN0 

XN0 

-2 

3xno 

3xno 

2xno 

2xno 

2xno 

kno 

X0+ 

X0+ 

x0+ 

0 

(3Xq+  -2) 

3X0+ 

2Xq+ 

2Xq+ 

2X0+ 

Ko+ 

d 

a* 

XN+ 

2nu 

Xn+ 

XN+ 

0 

3%+ 

(3V 

-2)  2Xjj+ 

2Xn+ 

2Xn+ 

* 

Kn+ 

XN0+ 

XN0+ 

XN0+ 

0 

3xno+ 

3%0* 

(2X  . 

NO 

-2)  2Xnq+ 

2xno+ 

kno'1 

x + 
°2 

x + 
°2 

X + 
°2 

a 

”<■£ 

2X  + 
°2 

(2X  . -2) 
°2 

2X  + 
°2 

X + 
N2 

x + 
N2 

X + 
N2 

0 

3X»5 

5X»; 

2X  + 
N2 

2X  + 
N2 

(2X  . -2] 
N2 

28 


NOLTR  74-182 


The  temperature  dependence  for  the  equilibrium  constants 
corresponding  to  these  reactions  are  presented  in  BORTNER  and  BAURER 
(1972)  . 

Some  brief  comments  concerning  the  selection  of  reaction  rates 
are  in  order.  The  reaction  rates  presented  in  BORTNER  (1969) 
represent  empirical  fits  made  to  accommodate  both  low  (T<500°K)  and 
high  (T>2500°K)  temperature  data  when  possible.  BORTNER  (1969)  and 
BORTNER  (1963)  reaction  rate  data  were  compiled  for  reentry  flow 
fields  and  are  appropriate  for  airshock  chemistry  calculations.  The 
corresponding  reaction  rates  presented  by  BORTNER  and  BAURER  (1972) 
are  concerned  with  low  temperature  data.  BORTNER  and  BAURER  (1972) ; 
however,  provide  a more  inclusive  list  and  present  the  most  recent 
estimates.  In  many  cases  the  reaction  rates  compiled  in  BORTNER  and 
BAURER  (1972)  represent  estimates  determined  by  a committee  of  workers 
knowledgeable  in  the  field. 

All  positive-ion  charge  exchange  reaction  studies  were  performed 
at  low  temperatures.  BORTNER  (1969)  suggests  an  uncertainty  factor 
of  (T/300)±1.5  for  which  the  exponent  value  1.5  is  an  estimate. 
Reaction  rates  for  positive-ion  charge  exchange  reactions  are  selected 
from  BORTNER  and  BAURER  (1972)  for  the  present  calculations  and  the 
temperature  factor (T/300 )~1 • 5 is  used  unless  otherwise  specified  by 
the  references.  For  dissociation,  rearrangement,  and  ionization 
reactions  the  reaction  rates  are  taken  from  BORTNER  (1969)  when 
available;  otherwise,  they  are  taken  from  BORTNER  (1963) . However, 
one  class  of  dissociation  (or  association)  reactions  are  taken  from 
BORTNER  and  BAURER  (1972) ; the  positive-ion  - neutral  dissociation 
reactions  are  not  considered  in  the  earlier  two  reaction  rate 
compilations . 

Each  reaction  rate  is  determined  on  an  individual  basis,  but 
generally  the  guidelines  discussed  above  are  used  to  select  the 
reaction  rates.  Parameters  from  BORTNER  (1963)  and  BORTNER  (1969) 
are  normalized  to  300°  K and  reaction  rates  defined  in  the  manner 
kRR  <lxlO-A  are  re-defined  as  kRR=lxlO-B  where  B=A+1. 
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TABLE  A-l 


REACTIONS  AND  REACTION  RATE  CONSTANTS 


Dissociation  Reactions 
Reaction 


0,  + 0„ 

2 2 


0 + 0 + 0, 


°2  + N2 


± 0 + O + N- 


°2  + ° 


->•  0 + 0 + 0 


LO 


o2  + X 


-*•  0 + 0 + X 


N„  + N„  N + N + N 

2 2 *■  2 


N„  + N -*■  N + N + N 
2 «- 


Reaction  Rate 

k = 1.2  x 10-7±’5 (T/300) -,83exp(-59400/T) 
kr  = 3.2  x 10-33i*  7 (T/300) -1*  3 
k = 5.0  x 10_7i*7 (T/300) -1exp(-59400/T) 

kf  = 8.34  x 10~8±-5 (T/300)-1*8exp(-59400/T) 

k = 1.15  x 10” 7i*5( T/300)-’ 75exp (-113200 /T) 
k = 9.6  x I0~6i*6 (T/300) ~1*5exp(-113200/T) 


Reference* 

A2 

A2 

A2 

A2 

A2 

A2 
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TABLE  A-l  CONTINUED 


Reaction 


Reaction  Rate 


Re ference* 


NT2  + X •*  N + N + X 


kf  = 4.47  x 10  8-' ' (T/300)  * 82exp (-113200/T)  A2 


NO  + N2  j N + 0 + N2 


kr  = 1.1  x 10  32±*5 


A2 


NO  + X -+  N + O + X 


kf  = 1.5  x 10  9-'  'exp  (-75400/T) 


A2 


n2+  + X 


-*•  NT  + N + X 


k^  = 1.0  x 10“29±2 (T/300) -1-5+1*5 


A3 


LO 


+ 


^ 0 + X 


+ 0 + 0 + X 


k = 1.0  x 10 
r 
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A3 


NO+  + X N+  + 0 + X 


kr  = 1.0  x 10“29i2 (T/300)-1*5-1'5 


A3 


NO+  + X ± 0+  + N + X 


kr  = 1.0  x 10 


-31 


A3 
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TABLE  A- 1 CONTINUED 


Rearrangement  Reactions 
Reaction 


0 + N2  £ N + NO 

0 + NO  - N + 0„ 

«-  2 


Reaction  Rate  Reference* 

kf  = 1.0  x 10_1°— * 3exp (-38000/T)  A2 

k = 3.7  x 10-13i*  3 (T/30  0)  1 ' 5exp (-19 10  0/T)  A2 


UJ 

u> 
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TABLE  A-l  CONTINUED 


Ionization  Reactions 
Reaction 


Reaction  Rate 


N + 0 


NO+  + e 


kr  = 4.62  x 10"7i*5 (T/300)-1*5 


N + e N + e + e 


= 4.87  x lO-30^1 (T/300)  1 


N + X ->■  N+e  + X 


, , .n-36+2 

k = 3.31  x 10  — 

r 


u> 

^ 0 + e 


t O + e + e 


kr  = 2.57  x 10"30!1 (T/300)-1 


O + X ■* 


O + e + X 


kr  = 1.9  3 x 10  36i2 


NO  + e NO  +e  + e 


kr  = 1.06  x 10  27i1(T/300)  1,5 


NO  + X -*■  NO  + e + X 


-34+2  - 5 

kr  = 7.96  x 10  J - (T/300)  ' 


N + N + N_"t’  + e 


k^  = 4.79  x 10-7!1 (T/300) -’5 


Refe  rence* 
A2 

A1 

Al 

A1 

Al 

Al 

Al 

Al 
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TABLE  A-l  CONTINUED 


Ionization  Reactions 
Reaction 

0 + 0 -+  0/  + e 

«-  2 


UJ 

Ln 


N2  + N 

N2  + e 

N2  + X 

0 + N 

2 

o2  + 0 
02  + e 
o2  + X 


2 


2 


2 


N2+  + e + N2 

N2+  + e + e 

N2+  + e + X 

0 + + e + N 
2 2 

02+  + e + 02 
02+  + e + e 

J 

02+  + e + X 


Reaction  Rate 
= 2.77  x 10-7— 1 (T/300) -1 

kr  = 5.2  x 10"35±2  (T/300)  ,5 

kr  = 1.77  x 10_29±1 (T/300)-1 

kr  = 3.58  x 10-38±2 (T/300)1*5 

k = 1.6  x 10~36±2 
r 

kr  = 3.2  x 10“36±2 

k^  = 1.17  x 10"29±1  (T/300)-1 

k = 1.37  x 10"35±2  (T/300)-*5 
r 


Reference* 

A1 

A2 

A2 

A1 

A2 

A2 

A2 

A1 
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TABLE  A-l  CONTINUED 


co 

T> 


Positive-Ion  Charge  Exchange 
Reaction 

N+  + O ± N + 0+ 


N+  + 02  1 NO+  + O 


N+  + NO  -*■  N + NO+ 


0+  + N0  NO+  + N 

Z Mr 


0+  + NO  -*■  0 + NO+ 


N2+  + N + N + N+ 


N + + 0 ->  N„  + 0+ 
2 -<-2 


N2+  + 0 + NO+  + N 


Reaction  Rate 
= 1.0  x 10"12  (T/300)i1,5 

= 3.0  x 10“10  (T/300)i1,5 

= (8  + 2.4)  x 10-10  (T/300)i1*5 

= 1.2  x 10-12 (T/300)  0 * 5 

= 1.3  x 10 ~13 (T/300) i1 ' 5 

= 1.0  x 10"12  (T/300)i1,5 

= 6.0  x 10-12  (T/300)i1,5 

= 2.5  x 10-10i*  7 (T/300)  — 1 *5 


Refe  re  nee* 
A3 

A3 

A3 

A3 

A3 

A3 

A3 

A2 
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TABLE  A-l  CONTINUED 


u> 

-j 


Positive-Ion  Charge  Exchange 
Reaction 


N2+  + °2  : N2  + °2+ 


N + + 0„  N0+  + NO 

2 2 «* 


N*  + NO  -*•  N_  + N0+ 
2 «-  2 


N+  + 02  j N + 02+ 


0+  + o2  + 0 + o2+ 


0+  + NO  -*■  0 + + N 

<-  2 


0 + + N -*■  NO+  + NO 
2 2 ■*- 


0 + + N -*■  NO+  + 0 


Reaction  Rate 

i i « . . — 10  + . 3 . . +1.5 

kf  = 1.0  x 10  — (T/300)  — 

kf  = 1.0  x 10“17i2  (T/300)-1'5 

k = 5.0  x 10"10i*3(T/300)i1'5 

kf  = (3.0  + 1.5)  x 10~10  (T/300)!1'5 

kf  = 2.0  x 10_11(  T/300)  ~ ’ 5 

k = 1.0  x 10-13 (T/300) -1 ' 5 

+1 

k = 1.0  x 10-16-3 (T/30 0) — 1 * 5 
kf  = (1.8  + .6)  x 10~10  (T/300)!1-5 


Reference* 

A2 

A2 

A2 

A3 

A3 

A3 

A3 

A3 
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TABLE  A-l  CONTINUED 


Positive-Ion  Charge  Exchange 
Reaction 


02+  + NO  j 02  + N0+ 


NO+  + N + X -*■  0+  + N-  + X 

2 


NO+  + 0 + X 0+  + NO  + X 


Reaction  Rate 

kf  = (6.3  + 2.4)  x 10"10  (T/300)i1,5 
kr  = 2.8  x 10"28(T/300)"1*5t1’5 
kr  = 1.0  x 10"29i2 (T/300)_1*5i1,5 


u> 

TO 


Reference* 

A3 

A3 

A3 


*Reference  Notation: 
BAURER  (1972). 


A1  - BORTNER  (1963) , A2  - BORTNER  (1968) f A3  - BORTNER  and 
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APPENDIX  B:  VIBRATIONAL  RELAXATION  EQUATION 


The  relaxation  equation  for  vibrational 
below* 


energy  is  presented 


dE 

v 

dx 


(Ev) 


- E 


VIB.  EQ. 


v 


Tp  (U  - U) 


k© 


v 


(N+l)k0 


v 


0 /T 

(e  v'  m_i ) 


a»ijVV.u 


- E 


v (X) 


+ (2  k0v  N - EV)  (X) 


Q- 


1_  /d(X)\ 

X)  \ dx  ) 


REC 


(XN  )nu 


DISS 


There  is  a separate  equation  for  each  species  which  undergoes 
vibrational  relaxation.  Vibrational  relaxation  is  considered  for 
two  species  in  this  model,  O2  and  ^ . The  other  molecular  species 
are  assumed  to  be  formed  in  vibrational-translational  equilibrium. 


The  first  term  in  the  above  equation  considers  collisional- 
vibrational  energy  exchange.  Values  for  characteristic  relaxation 
times,  t,  are  taken  from  MILLIKAN  and  WHITE  (1963).  The  simple 
harmonic  oscillator  model  is  used  for  computing  the  vibrational 
energy  as  a function  of  either  the  local  translational  temperature 
(EV)VIB  eq  (T)  or  the  local  vibrational  temperature  Ev  (Tv) . 

The  middle  two  terms  correct  the  vibrational  relaxation  rate  by 
including  the  effects  of  the  dissociation  process  (TREANOR  and 
MARRONE  (1962)).  The  more  energetic  vibrators  are  more  easily 
dissociated.  This  process  lowers  the  average  vibrational  energy  in 
the  molecules  that  remain.  However,  the  assumption  that  dissociation 
does  not  perturb  the  Boltzmann  distribution  is  maintained.  An  addi- 
tional correction  is  included  to  take  into  account  the  change  in  the 
number  of  oscillators  during  the  dissociation  process . 


The  last  term  in  the  relaxation  equation  considers  electron- 
vibrational  energy  exchange  which  only  applies  to  N2  and  will  be 
discussed  in  connection  with  the  electron  energy  equation. 

*See  6.  LIST  OF  SYMBOLS  of  text. 
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APPENDIX  C:  ELECTRON  ENERGY  EQUATION 

The  form  of  the  electron  energy  equation  used  is* 


dTe 
“ 3x 


5 

3 T 

rn  uk 
2 e 


T dnQ 
e e 

n dx 
e 


with  the  following  electron  energy  flux  terms. 


Q1  = Z (6. 85 3x10 ~4 3 ) (T-T|/2  (£n ( (7.675xl07)  + 1 .0) (X) ± (X  )n2 . 

i m.  T ' e 

i e 

is  the  energy  gained  by  electrons  through  elastic  collisions 
with  ions  (BIBERMAN  and  YAKUBOV  (1965) ) . Summation  is  over  all  ions 
i in  the  air  chemistry  model.  Electron-neutral  elastic  collision 
energy  exchange  is  negligible  in  comparison  and  is  not  included. 

q2  “ -5Eii<Dea>in2 


Q2  is  the  energy  lost  by  electrons  through  electron-neutral 
collisional  ionization  reactions.  (Dea^i  t*le  net  production  rate 
(CM^/SEC)  , forward  minus  reverse,  for  the  i^  electron-neutral 
collisional  ionization  reaction  and  Ej.  is  the  corresponding  ioniza- 
tion energy;  i is  summed  over  the  electron-neutral  collisional 
ionization  reactions  in  the  air  chemistry  model. 


,-9 


-10  /T 


Q3  = (4.5x10  )e 

• (Xe)(XN2)n2, 


k0v  (e  (TVN2)  -e  (Tp)  ) (e(TVN2)+e(T  )+l) 
N2 


where 


<®v>n/T  Y1 


e(T)  = \e  -lj 


*See  6.  LIST  OF  SYMBOLS  of  text. 
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Q3  is  the  energy  gained  by  electrons  through  collisions  with  N2 
molecular  oscillators  (BIBERMAN,  MNATSAKANYAN,  and  YAKUBOV  (1971)). 

In  general,  the  direct  excitation  of  vibrational  levels  of  diatomic 
molecules  by  collision  with  slow  electrons  is  not  efficient  because 
of  the  differences  in  masses  of  the  collision  partners.  However,  for 
N2  the  slow  electron  is  captured  forming  an  unstable  negative  ion. 
Following  the  decay,  the  final  vibrational  quantum  number  may  be 
several  steps  above  the  initial  quantum  number  which  effects  a con- 
siderable change  in  the  electron  energy.  This  mechanism  also  applies 
to  CO  which  is  not  included  in  the  present  air-chemistry  model.  Q3 
is  in  the  N2  vibrational  relaxation  equation. 


54=IkrV£n2-rkTeVr": 


Q4  is  the  energy  gained  by  electrons  from  associative  ionization 
reactions  (ZHELEZNYAK  and  MNATSAKANYAN  (1968)).  (kaj)f  and  (kAI)r 
are  the  net  forward  and  reverse  reaction  rates  (CM^/SEC)  for  associa- 
tive ionization  summed  over  reactions  in  the  air  chemistry  model. 

Q5  = -n2  (Xe)  Pe(Te)  (E0X0  + EnXn  + E0  X0  + E0+X0  + EN+XN+ 


+ E ,X  , + E .X  , ) 

0+  0+  n+  nJ 


t 


where  Pe(Te)  = a£ve  = 6.212  x 10“^^  Te^/^.* 


Q5  is  the  energy  lost  by  electrons  maintaining  the  quasistationary 
population  of  excited  states.  The  electron  excitation  energies  for 
species  Nj,  NO,  and  NO+  are  negligible  in  comparison  with  contributions 
from  the  other  species . 


*An  approximate  value  of  ae  = 10  ^CM2  was  selected  (ZEL'DOVICH  and 
RAIZER  (1966)). 
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APPENDIX  D:  ENTHALPY  EQUATION 


Enthalpy  is  defined  as* 


h(T,  T , T , X.)  = |(1-X  )RT  + |x  R T 
v e i 2 e 2 e e 


+ 


R(0V)  . 
1 


j/T. 


+ RT 


e 


E 


i 


i 


The  summations  are  given  to  be : j - molecules  and  molecular  ions , 

% - all  chemical  species,  i - electronic  levels.  The  vibrational 
model  used  is  the  harmonic  oscillator.  Separate  vibrational  tempera- 
tures are  considered  for  Oo  and  N2  vibrators  whereas  all  other  species 
are  evaluated  at  the  translational  temperature.  All  rotational  energy 
modes  are  assumed  to  be  excited  at  the  classical  limit  in  equilibrium 
with  translation,  and  the  distribution  of  electronic  states  corresponds 
to  the  free  electron  temperature;  that  is,  the  bound  and  free  electrons 
are  in  equilibrium.  The  thermodynamic  constants  are  taken  from  BORTNER 
and  BAURER  (1972)  . 


*See  6.  LIST  OF  SYMBOLS  of  text. 
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n/(n)|viAX 

V^e^MA 

p/(P>MAX 


FIG.  1.  TOTAL  NUMBER  DENSITY,  ELECTRON  NUMBER  DENSITY,  AND  PRESSURE  PROFILES  IN  RELAXATION 
REGION  FOR  M - 16.47,  p = 5 TORR,  T = 300° K 
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n/*n>MAX 


1.0 


0.8 


0.6 


(n)M AX  = 2.4x10'  8cm“ 


LABORATORY  TIME  BEHIND  SHOCK  FRONT-MICROSECONDS 


FIG.  2.  NEUTRAL  SPECIES(02,  N2,  NO,  O,  N)  PROFILES  IN  RELAXATION  REGION  FOR  M = 16.47, 
p = 5TORR,T  = 300°K 
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FIG.  3.  CHARGED  SPECIES  (N2t  N0+,  0+,  N+)  PROFILES  IN  RELAXATION  REGION  FOR 
M = 16.47,  p = 5 TORR,  T = 300° K 


t/(T)max 


FIG.  4.  TEMPERATURE  (TRANSLATION,  ELECTRON,  O2  VIBRATION,  N2  VIBRATION)  PROFILES  IN 
RELAXATION  REGION  FOR  M = 16.47,  p = 5 TORR,  T = 300°K 
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FIG.  5.  TOTAL  NUMBER  DENSITY  AND  ELECTRON  NUMBER  DENSITY  PROFILES  IN  RELAXATION 
REGION  FOR  M = 26.76.  p = 5 TORR,  T = 300° K 
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LABORATORY  TIME  BEHIND  SHOCK  FRONT-MICROSECONDS 


FIG.  6.  CHARGED  SPECIES  (02+,  N2+,  NO+,  0+,  N+)  PROFILES  IN  RELAXATION  REGION  FOR 
M = 26.76,  p - 5 TORR,  T = 300°K 


NOLTR  74-182 


P0TL  (TORR/iSEC) 


10"  F 


103  F 


10;  - 


10  r 


1 r 


10  1 z 


10J  r 


10  3 L 
6 


EXPERIMENTAL  DATA 
o SHAFER  AND  FROHN  (1972) 

• LIN.  NEAL  . AND  FYFE  (1962)  TAKEN  FROM  BIBERMAN,  MNATSAKANYAN, 


AND  YAKUBOV  (1971) 

• FROHN  AND  DE  BOER  (1969)  TAKEN  FROM  STRICKER  AND  LOW  (1972) 
a MANHEIMER-TIMNAT  AND  LOW  (1959) 
a STRICKER  AND  LOW  (1972) 

BASHENOVA  AND  LOBASTOV  (1972) 

NIBLETT  AND  BLACKMAN  (1958) 

WILSON  (1966) 

AIRSHOCK  RELAXATION  COMPUTER  CODE  CALCULATION 


J l I l I I I I I I I I I I I I I 

8 10  12  14  16  18  20  22  24  26  28  30  32  34  36  38  40 

MACH  NUMBER 


FIG.  7 DEPENDENCE  ON  MACH  NUMBER  OF  THE  PRODUCT  OF  IONIZATION  RELAXATION 
TIME  IN  LABORATORY  COORDINATES  AND  UPSTREAM  AMBIENT  PRESSURE  IN  AIR 
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FIG.  8.  DEPENDENCE  ON  MACH  NUMBER  OF  02  HALF-DISSOCIATION  LENGTH  IN  AIR 
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AIRSHOCK  RELAXATION 
COMPUTER  CODE 
CALCULATION 

EMPIRICAL  FIT 
TO  EXPERIMENTAL 
DATA-MILLIKAN 
AND  WHITE  (1963) 


I I I I I I I 

1 0.02  0.04  0.06  0.08  0.10  0.12 

(T°Kr,/3 

FIG.  9.  DEPENDENCE  ON  TEMPERATURE  OF  THE  PRODUCT  OF  VIBRATIONAL  RELAXATION 
TIME  IN  PARTICLE  COORDINATES  AND  DOWNSTREAM  PRESSURE  FOR  02  AND  N2 
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